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Abstract 

We use a mapping of the multiband Hubbard model for CuOs chains in 
i?Ba2Cu306+a; {R =Y or a rare earth) onto at — J model and the description 
of the charge dynamics of the latter in terms of a spinless model, to study 
the electronic structure of the chains. We briefly review results for the optical 
conductivity and we calculate the quantum phase diagram of quarter filled 
chains including Coulomb repulsion up to that between next-nearest-neighbor 
Cu atoms V2, using the resulting effective Hamiltonian, mapped onto an XX Z 
chain, and the method of crossing of excitation spectra. The method gives 
accurate results for the boundaries of the metallic phase in this case. The 
inclusion of V2 greatly enhances the region of metallic behavior of the chains. 
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There is consensus in that the electronic structure of it!Ba2Cu306+j; {R —Y or a rare 
earth) can be separated into that of the two Cu02 planes per unit cell which become su- 
perconducting under doping, and that of the Cu02+a; subsystem, in which CuOs chains are 
formed for oxygen content x > 0.5 and low temperatures. [1] The electronic structure of the 
CuOa chains is crucial because it controls the doping of the superconducting Cu02 planes. 
The dependence of the superconducting critical temperature Tc with annealing, [2] combined 
with Raman measurements [3] and persistent photoconductivity experiments [4,5] show an 
intimate relation between the oxygen ordering in the CuOa; planes and Tc [6] : oxygen order- 
ing along chains increases the amount of two-fold and four-fold coordinated Cu atoms at the 
expense of three-fold coordinated ones, and leads to an increase in the hole doping of the 
superconducting Cu02 planes. Detailed calculations of the relation between electronic and 
atomic structure in i?Ba2Cu306+a;, together with a simple explanation of the above facts 
valid in the strong coupling limit were presented. [1] These results show the relevance of 
interatomic Coulomb interactions. In addition even near the optimum doping (~ 1/5 holes 
per Cu atom in the planes), the average distance between carriers is of the order of two 
lattice parameters of the planes suggesting that interatomic repulsion at smaller distances 
are screened only partially. 

Several pieces of evidence suggest that the CuOa chains are insulating. For example 
PrBa2Cu307 is semiconducting, [7] and the contribution of the CuOa chains to the optical 
conductivity cr(a;) is very similar in this compound [8] and in superconducting YBa2Cu306+a;, 
[9,10] displaying a broad peak near u ~ 0.2 eV and a slowly falling tail at higher frequen- 
cies. Also, charge modulations observed by scanning tunneling microscopy (STM) were 
interpreted in terms of a charge density wave and a gap in the spectrum of the chains. [11]. 
Finally in the explanation of Fehrenbacher and Rice of the suppression of superconductiv- 
ity upon substituting Y by Pr in YBa2Cu307, they propose that the holes which dope the 
superconducting Cu02 planes in YBa2Cu307 are displaced towards a hybrid Pr-0 state in 
PrBa2Cu307. [12] This implies a shift in the Fermi level of about 0.25 eV according to their 
parameters, while the authors assume that the hole occupation of the CuOs chains is 0.5 in 
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both cases. Thus, this explanation seems to require a gap in the chains to be consistent. 

However, all he above data can also be consistently explained assuming intrinsically 
metallic chains cut by ~ 5% of defects or oxygen vacancies {x ~ 0.95), which is usual in these 
systems. [13] The appropriate multiband model for CuOa chains was mapped numerically 
into a t — J model with t ~ 0.85 eV and J~0.2 eV. With these parameters, the decrease 
in the occupation of the chains upon replacing Y by Pr is only 0.05. Taking into account 
that the charge dynamics of the model can be described up to a few percent by a spinless 
model even for J/t — 0.4, [13,14] charge modulations and the optical conductivity can be 
explained. [13] In particular, the lower energy part of the latter is given by: 

AB 

with 

A = — 2t7rln(l — c) sin kp; B = —{1 — c^)tsin kp, 

and c = 1 — a; is the concentration of oxygen vacancies. The experimental results were 
fitted using A — 0.35 eV. [13] The resulting optical conductivity below 0.4 eV is shown 
in Fig. 1. Note that in spite of the metaUic character of the chains, as a consequence of 
the oxygen defects, cr(a;) has a pseudogap at low energies. The experiments cannot confirm 
this due to large errors for u; < 0.1 eV. [8-10] However, recent STM studies of the local 
density of states detect a pseudogap of about 25 meV and numerous intragap resonances. 
[15] The latter might be explained by the effect of defects on superconductivity in the chains 
induced by proximity, [16] but also in principle by eigenstates of long finite metallic chains. 
Unfortunately, the local density of states of the one-dimensional t — J model depends also 
on the spin wave function and cannot be described solely by spinless fermions. [17,18] 

The natural candidate to open a gap in the effective t — J model for the CuOa chains 
is the nearest-neighbor repulsion Vi. Keeping the assumption that the charge dynamics is 
described by a spinless model, one expects that a gap opens for Vi > 2t ^ 1.7 eV. [19,20] 
If the Coulomb repulsions were completely unscreened Vi ~ e^/b = 3.6 eV, where b is 
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the lattice parameter along the chains. Recently Seo and Ogata showed that inclusion of 
next-nearest-neighbor repulsion V2 enhances the range of stability of the metallic phase, 
calculating the gap as a function of V2. [21] 

We calculate the phase diagram of the spinless model, including Vi and V2 using the 
method of crossing of excitation levels. [22-24] Actually, the mapping of the energy of the 
one-dimensional t — J model into that of a spinless model is strictly valid only for J = 0, 
[17,20,21] but we expect it to be a very good approximation for J/t < 0.4. [13,14] The 
advantage of the method of level crossings, briefly explained below, over previous approaches 
[21,25] is the accuracy that can be achieved for the phase boundaries. This has been shown 
for example in its application to the Hubbard model with correlated hopping [26,27] in 
comparison with exact results. [28] 

In standard notation, the model is: 

H = ^[t(4+iCi + H.c.) + ViUiUi+i + V2nini+2\, (2) 

i 

with Tij = c'-Cj. Using a Jordan- Wigner transformation = exp(i7r Ylii<j — i^j~V^ 

Sj — rij — 1, the model can be mapped into an XXZ model with next- nearest- neighbor 
antiferromagnetic Ising interaction: 

H = J^^MSfS:^, + SfSl,) + A.StSt^, + A^SfSt^,], (3) 

i 

where is the /S component of the spin-1/2 operator at site i, Ji — 2t and Aj — Vj. [20,21] 
A successful approach to describe the qualitative properties of one-dimensional strongly 
correlated systems is bosonization followed by a renormalization group procedure. This 
procedure usually terminates at a fixed point, which determines the properties of the system 
for the initial parameters given. A phase transition occurs when the flow goes towards a 
different fixed point. Since the renormalization group is a weak coupling approach, the 
phase boundaries are not given accurately by the method for large interactions. The basic 
idea of the method of level crossings is to combine numerical calculations of excitation levels 
with basic knowledge on the properties of these fixed points. The more interesting phase 

4 



transitions involve one fixed point which is scale invariant. This is for example the case of 
the XX Z model with next-nearest-neighbor interactions studied by Nomura and Okamoto. 
[22] The spin fluid phase of Eq. (3) (which corresponds to the metallic phase of Eq. (2)), like 
that of an ordinary Heisenberg model is characterized by a scale invariant fixed point. [22] 
Then, using conformal field theory one can relate the excitation energy which corresponds 
to some operator Ai at site i (for example a spin flip , S^), to the dependence of the 
correlation functions of this operator with distance d, for large d: 

E4L)-E,{L)^^-^, {A+,A,)^J^. (4) 

Here L is the length of the system, v the spin- wave velocity, Eg{L) the ground state energy, 
Ea{L) the lowest energy in the adequate symmetry sector (connected to the ground state 
by Ai) and xa the critical dimension for the excitation A. Since the dominant correlations 
at large distances determine the nature of the thermodynamic phase, a phase transition is 
determined by the crossing of excited levels for different symmetry sectors. 

In the present problem, the relevant quantum numbers which determine the symmetry 
sector are total wave vector K, total spin projection S^, parity under inversion P and parity 
under time reversal T . We have restricted our calculations to number of sites L multiple 
of four to avoid frustration of the phase which we call AFII (see below). For these sizes, 
the quantum numbers of the ground state are always the same in the region of parameters 
studied. They are listed in Table I, together with the quantum numbers of the first excited 
state of each phase. Our main interest are the boundaries of the spin fluid phase of the spin 
model Eq. (3) which corresponds to the metallic phase of Eq. (2). With increasing Ai (A2) 
there is a continuous transition to an insulating Neel ordered (dimerized) phase. [22] The 
Neel ordered phase, which we call antiferromagnetic I (API) for maximum order parameter 
has a spin ordering titi ■■■ corresponds to a charge ordering 1010... in the original model 
Eq. (2). The dimer phase has a gap which is exponentially small near the metallic phase [22]. 
This renders it very difficult to detect the transition with alternative numerical methods. 
[27] The transitions between any two of these three phases were determined accurately from 
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the corresponding crossing of excited levels (see Table I). In addition, with increasing A2, 
we expect a transition from the dimer phase to an AFII a phase with long range order 
ttii ••• (corresponding to charge ordering 1100...). This transition cannot be detected by 
crossing of first excited states. Since it involves two insulating phases, it is not described 
by a scale invariant theory and is also beyond our scope. For the sake of completeness we 
have drawn a tentative dimer- AFII boundary using the rough criterium that the system is in 
the AFII phase when the ground state correlation function (calculated deriving the energy 
using Hellmann-Feynman theorem) {Sf Sf j^2) < —1/8. For the other transitions, we have 
calculated the transition points in systems with L = 12, 16 and 20 sites. According to field 
theory predictions for large enough L, these points plotted as a function of 1/L^ should he 
on a straight line. [22,29] We have verified that this is the case for the three transitions with 
high accuracy. The linear fit provided the transition point extrapolated to 1/L^ — > 0, and 
is error. The error is below 1% in all cases, confirming the validity of the method in the 
present case. 

The resulting phase diagram is shown in Fig. 2. For V2 = A2 = 0, the known exact 
results [19,20] are reproduced: there is a transition from the spin fluid (metallic) phase to 
the AFI (charge density wave) phase at Ai — Ji (Vi — 2t). Another known limit is the 
classical one Ji — (t — * 0), for which there is a transition between both AF phases at 
A2 = Ai/2 (V2 — ^i/2). Our results are consistent with this hmit. However, there is a strip 
of width ~ Ji = 2i of a dimer phase between both AF (charge ordered insulating) phases. 
This is reminiscent of the physics of the ionic Hubbard model, for which a strip of a dimer 
phase of width ~ 0.6t in the strong coupling limit, separates the band insulating and the 
Mott insulating phases [24] dye to the charge fluctuations that still remain in the strong 
coupling limit. 

In qualitative agreement with previous calculations, [21] we obtain that the addition 
of V2 greatly enhances the range of stability of the metallic phase of the CuOs chains in 
i?Ba2Cu306+2:. For unscreened interactions Vi = 2V2 — 3.6 eV. Using Ji = 2t = 1.7 eV, 
one can see from the phase diagram, that the system falls in the metallic phase even in 
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this extreme case. Instead, if V2 were neglected the chains would be in an insulating charge 
ordered state for the same t and Vi. 

A numerical mapping of the appropriate multiband Hubbard model for the chains to a t — 
J model indicates that J/t < 1/4 [13]. It is reasonable to expect that turning J to zero does 
not change sustantially the phase diagram. For J — 0, the mapping to the spinless model is 
exact and our results lead to the conclusion that the CuOa chains in i?Ba2Cu306+a; (-R =Y 
or a rare earth) are intrinsically metallic. Observed charge modulations are likely due to 
Priedel oscillations induced by defects, hke O vacancies. These defects or superconductivity 
induced by the CUO2 planes can also lead to the observed pseudogap behavior. 
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grateful to the National Research Council (CNPq), Brasil, by their financial support. A.A. 
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TABLE I 





K 




p 


T 


ground state 








1 


1 


exc. spin fluid 


TT 


±1 


-1 




exc. AFI 


TT 





-1 


-1 


exc. dimer, AFll 


TT 





1 


1 



Quantum numbers of the ground state and the first excited state of the different phases 
for L multiple of four. 



FIGURE CAPTIONS 

Fig.l. Low energy part of the optical conductivity of CuOa chains for A — 0.35 eV 
(see Eq. (1)). 

Fig.2. Phase diagram of the effective model for CuOa chains (Eq. (2) or (3)) as a 
function of Ai/Ji = Vi/2t and A2/J1 = ¥2/21. 
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